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Time-dependent density-functional theory (TDDFT) often successfully reproduces excitation en¬ 
ergies of finite systems, already in the adiabatic local-density approximation (ALDA). Here we show 
for prototypical molecular materials, i.e., oligothiophenes, that ALDA largely fails and explain why 
this is so. By comparing TDDFT with an in-depth analysis based on many-body perturbation the¬ 
ory, we demonstrate that correlation effects crucially impact energies and character of the optical 


excitations not only for molecules of increasing 
for isolated small molecules. We argue that only 
correlation effects, can reproduce optical spectra 
gas phase to crystal structures. 

PACS numbers: 71.15.Mb, 71.35.Cc, 78.40.Me 

I. INTRODUCTION 

Optical excitations in organic materials are strongly 
dominated by many-body effects. Electron-electron 
interactions determine the electronic structure, and 
electron-hole (e-/i) correlations rule the excitation pro¬ 
cess. A methodology that is able to consistently cap¬ 
ture the features of molecular materials, from single 
molecules to their condensed phases, including organic 
crystalsad sorba te systems]^ hybrid materials and 
nanostructures, P^° l ^^ l is an essential prerequisite to pre¬ 
dict their excited-state properties. Many-body pertur¬ 
bation theory (MBPT) represents the state-of-the-art 
method to calculate optical excitations in solids .1^^ The 
GW approaclf^ gives quasi-particle (QP) energies and 
the solution of the Bethe-Sapleter equation (BSE)pIl3l 
yields excitation energies and wavefunctions of the e- 
h pairs. Although, in principle, GIU-I-BSE can be ap¬ 
plied to any material, it is computationally too demand¬ 
ing for many systems of technological interest. Quantum 
chemistry offers powerful tools, such as Coupled Clusters 
(CC) and Configuration Interaction Singles and Doubles 
(CISD),Ii^ to accurately compute optical properties of 
molecules. Again, only small systems can be treated with 
these methods, due to their huge computational costs. 

Since the turn of the century, time-dependent density- 
functional theory (TDDFT) has gained increasing 
popularity,li^ due to its remarkable ability in reproduc¬ 
ing optical spectra of small molecules and clusters, with 
relatively low computational effort This success is 
mainly ascribed to the dominance of the long-range 
part of the Coulomb potential v over the exchange- 
correlation kernel {fxc)- v blue-shifts the absorption 
onset above the Kohn-Sham (KS) gap, and redistributes 
the oscillator strength (OS) to higher energies, provid¬ 
ing a good approximation for the optical spectra.^ How¬ 
ever, TDDFT suffers from severe drawbacks when deal¬ 
ing with extended systems. The spurious long-range 


length and in crystalline environment, but even 
high-level methodologies, which explicitly include 
of molecular materials with equal accuracy from 


behavior of stan dard fxc inhibits reproducing bound 
excitons in solids.l^^^ For the same reason, TDDFT 
is unable to describe charge-transfer-like excitations in 
molecular complexesTo overcome these limita¬ 
tions, new kernels have been developed including many- 
body effect ^21H351 exhibiting the correct long-range 

behavior.l2fiH^Unfortunately, these improvements mainly 
concern specihc classes of materials and/or excitations, 
and therefore do not often extend, in practice, the range 
of applicability of TDDFT. 

The goal of this work is to understand the role of 
many-body effects in the optical excitations of molecu¬ 
lar materials, from the gas-phase to crystals, and to clar¬ 
ify when and why TDDFT can be trusted. We adopt 
the adiabatic local-density approximation (ALDA),l^as 
the simplest and most common kernel of TDDFT. While 
more soph istica ted kernels can quantitatively improve 
the results, 1^^*^211111 the physical picture is already clear 
from ALDA.I^ We compare TDDFT with MBPT, with 
the aim to identify the most relevant contributions to 
the optical excitations and to understand to which extent 
TDDFT is able to reproduce them. To do so, we choose 
the family of oligothiophenes as prototypical example. 
Since only a few intense peaks characterize their UV- 
visible spectra,*^ they are ideally suited for this study. 
We investigate thiophene oligomers with an even num¬ 
ber of rings, up to 6, as well as the single thiophene ring, 
going from the isolated molecule to the crystal. 

II. THEORETICAL BACKGROUND AND 
COMPUTATIONAL DETAILS 

Optical absorption spectra presented in this work 
are computed from frequency-dependent linear-response 
TDDFT, as well as from MBPT, in the framework of 
GoWo-bBSE.EmiZl 

In TDDFT, optical excitations are calculated from the 
solution of the Dyson-like equation for the dynamic po- 
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larizability x = Xo + Xoi.v + fxc)x, where xo is the KS 
response func tion. In ALDA, the kernel is lo¬ 

cal and static.l232SI Por f^c = 0, X is computed within 
the random-phase approximation (RPA), while neglect¬ 
ing also V leads to the independent-particle approxima¬ 
tion (IPA), with X = Xo- The dielectric tensor e is calcu¬ 
lated as eG,G'(q,w) = (5 g,g' - i^G(q)XG,G'(q, w). 

Optical excitations are computed from MBPT through 
a two-step procedure.!^ First, the QP correction is ob¬ 
tained from GoIFo- Then, exciton energies are given by 
the solution of the BSE, in the matrix form: = 

E^A^. The effective two-particle Hamiltonian = 

H^'^ag is composed of three terms, which 

can be switched on and off, depending on the values of the 
coefficients Xx and Xc- For 7a;=7c=l (7a;=0, 7c=l) singlet 
(triplet) excitations are calculated. The diagonal term 
jjdiag accounts for single-particle transitions. The re¬ 
pulsive e-h exchange term , including the short-range 
Coulomb interaction, describes local-field effects (LFE). 
The direet term contains the screened Coulomb in¬ 
teraction, which determines the attractive e-h interac¬ 
tion. The imaginary part of the macroscopic dielectric 
function is related to the BSE eigenvectors A^, which 
correspond to the amplitude of the excitons. They also 
carry information about the character of the excitons and 
the weight of the single-particle transitions contributing 
to them. 

All calculations are performed with the exciting 
code,^ a computer package for density-functional theory 
and MBPT, implementing the all-electron full-potential 
augmented planewave method. The KS electronic struc¬ 
ture is computed using the Perdew-Wang local-density 
approximation (LDA) as the exchange-correlation (xc) 
functional.!^ In exciting, TDDFT and BSE are treated 
on the same footing,!^ enabling a direct comparison be¬ 
tween the results. Isolated molecules are considered in 
orthorhombic supercells, including at least 7 A of vacuum 
in each lattice direction. The internal coordinates are re¬ 
laxed, with a threshold of 0.025 eV/A for the interatomic 
forces. All the resulting geometries are flat. For the cal¬ 
culation of Xj 100 empty states are included. In GqWq 
calculations, the dynamically screened Coulomb poten¬ 
tial Wo is computed within RPA, including 100 empty 
states for the crystalline structures and 500 for the iso¬ 
lated molecules. For the solution of the BSE, 500 empty 
states are considered to compute the screened Coulomb 
interaction. The adopted computational parameters en¬ 
sure accuracy of 0.05 eV for the lowest-energy excitations 
in the spectra. 


III. RESULTS AND DISCUSSION 

A. Molecules of increasing length 

In Fig. [T] we show the optical spectra of bithiophene 
(2T), quarterthiophene (4T) and sexithiophene (6T). For 
the smallest oligomer, 2T, the ALDA spectrum is in ex- 


. 2T BSE 1 

“o 1 

— BSE 
— ALDA [ 
1 

LDA . GoWo] 

xm 

;4T 1 

A 

n i-A. , . -^ 

: 6T k 

III 

/i\ 

y ' \ 




0 1 2 3 4 5 6 7 

Energy [eV] 


FIG. 1: (Color online) Imaginary part of the macroscopic di¬ 
electric function of 2T, 4T, and 6T for light polarization par¬ 
allel to the long molecular axis (x). Results from ALDA, BSE, 
and an approximation of the e-h interaction term (Hq^^, see 
text) are presented. LDA and GolFo gaps are indicated by 
dashed lines. The molecules are shown as insets with C (S) 
atoms depicted in green (yellow) and H in black. A Lorentzian 
broadening of 0.1 eV is applied to all the spectra. 


cellent agreement with the GoWo+BSE result. The UV- 
visible region is dominated by a strong peak, polarized 
along the long (x) axis of the molecule,!^ and governed by 
the TT-TT* transition between the highest-occupied molec¬ 
ular orbital, HOMO (H), and the lowest-unoccupied 
molecular orbital, LUMO (L). Our Ending, indicating the 
lowest excitation energy of 2T at 4.05 - 4.09 eV (see Ta¬ 
ble]^, is in agreement wit h exp eriments,!^^!!^ and with 
quantum-chemistry results,which evaluate the peak 
between 4.05 and 4.11 eV. As the size of the oligomer 
increases, discrepancies between the GoHd+BSE and 
ALDA spectra emerge in both the energy of the first 
peak and the overall spectral shape. 

To understand the source of this disagreement, we in¬ 
spect in detail the results from MBPT concerning the 
first excitonic peak. This analysis is summarized in Ta¬ 
ble]^ Two types of many-body effects come into play: the 
QP correction to the electronic structure (S — 14c) and 
the e-h interaction. Like the fundamental gap computed 
from LDA and GqWo (Eff^°), also S - 14c, 

which corresponds to their difference, decreases with in¬ 
creasing oligomer length. Its effect is to blue-shift the 
absorption onset in absence of e-h correlation and LFE. 
From the solution of the BSE, the exciton binding ener¬ 
gies (E}f) are computed; for simplicity, we define E^ a s the 
difference between the excitation energy and E^f^° HeEl] 
Due to its bound and localized character, the lowest ex¬ 
citon in 2T has a large binding energy of almost 3 eV 
(Table |T|. The peak position arises from a partial can¬ 
cellation between E — 14c and E^. We label this dif¬ 
ference _ j-y; _ which is usually a 

positive quantity. This compensation should be mim- 
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jpLDA 

-^gap 

pGoWo 

^gap 

E-F,, 

Ej^BSE E-triplet 

Esingiet 

^MBPT 

^ALDA 

^TDDFT 

2T 

2.91 

6.82 

3.91 

4.65 2.19 

(-2.17) (-4.63) 

4.09 

(-2.73) 

1.18 

4.05 

1.14 

4T 

1.91 

4.21 

2.30 

2.70 0.73 

(-1.51) (-3.48) 

2.37 

(-1.84) 

0.46 

2.70 

0.79 

6T 

1.66 

3.21 

1.54 

2.04 0.28 

(-1.17) (-2.92) 

1.66 

(-1.54) 

0.01 

2.22 

0.56 

IT crystal 

4.42 

8.34 

3.92 

7.17 3.75 

(-1.17) (-4.59) 

5.80 

(-2.54) 

1.38 

5.44 

1.02 

IT expanded crystal 

4.63 

9.47 

4.84 

7.69 3.81 

(-1.78) (-5.66) 

5.77 

(-3.70) 

1.14 

5.63 

1.00 

IT molecule 

4.59 

9.85 

5.26 

7.14 3.78 

(-2.71) (-6.07) 

5.85 

(-4.00) 

1.26 

6.66 

2.08 


TABLE I: Energies relevant for the analysis of excitations in all systems considered in this work: Fundamental gaps, obtained 
by LDA and GqWo and their difference (E — 14c); exciton energies of the first triplet (Etripiet) and singlet {Esingiet) peak as 
well as from an approximation to the e-h interaction {Ejjbse, see text); = (E — 14c) + E), for the first bright exciton 

from BSE; _ j^alda _ with being the first peak energy from TDDFT. Exciton binding energies {Eb) 

are given in parenthesis. All energies are expressed in eV. 


icked by the ALDA kernel, in order to correctly yield 
the absorption features. For this purpose, we dehne 
^TDDFT^^ALDA _ , where denotes the 

lowest-energy bright excitation from ALDA. The closer 
^TDDFT jg ^MBPT ^ better the partial compensa¬ 
tion between E — 14c and Eb is reproduced by ALDA. In 
2T, J\TDDFT ^ ^MBPT. a_LDA mimics almost perfectly 
the blue-shift caused by the QP correction and the red- 
shift due to exciton binding. As the size of the oligomer 
increases, the spread between /\tddft ^mbpt 
comes larger. 

We note that we have to deal with two different issues 
here. One is a shortcoming of ALDA, the main focus of 
this work, which we discuss below. At this point, how¬ 
ever, we also need to critically assess the MBPT results in 
view of the GW starting-point problem. GolFo-l-BSE is 
known to underestimate excitation energies of gas-phase 
molecules^nd this shortcoming becomes more se¬ 
rious for large compounds. Indeed, we find the absorp¬ 
tion onsets of 4T and 6T to be red-shifted compared to 
the experi mental ones, assigned at 3.2 eV and 2.9 eV, 
respectively! ^^ * ^^ ! Th is issue is predominantly ascribed to 
the GqWo stepj ^° l ^^ l whose results may crucially depend 
on the underlying xc functionaland becomes very 
critical in case of large oligomers. In Ref. ESI GokFo gaps, 
computed on top of different xc functionals, are presented 
for a number of compounds, including oligoacenes. Over¬ 
all, LDA as starting point leads to a considerable under¬ 
estimation of the GokFo gaps compared to experiments, 
with the smallest oligomer described best. A system¬ 
atic analysis on the starting-point dependence of MBPT 
spectra in molecular materials from gas to crystal phase 
is definitely due, but goes beyond the scope of this work. 
This issue, although affecting the agreement with exper¬ 
iments, does not alter exciton binding energies and exci¬ 
tation character. As such, it does not impact the analysis 
that follows. 

We focus on the BSE results, which help us to clar- 



2T 

4T 

6T 

Eq 

H L (100%) 

H L (100%) 

H L (100%) 

triplet 

H ^ L (96%) 

H ^ L (89%) 

H ^ L (59%) 

H-1 ^ L-kl (27%) 

singlet 

H L (84%) 

H L (89%) 

H L (74%) 

H-1 L-kl (20%) 


TABLE II: Single-particle transitions contributing to the first 
peak of 2T, 4T and 6T. triplet and singlet excitation 

energies are shown. 


ify the shortcomings of ALDA. To this extent, we ana¬ 
lyze the nature of the low-energy excitations in terms of 
single-particle contributions, namely, whether they stem 
from one transition or from a combination of two or more. 
For this purpose, we indicate in Table the composi¬ 
tion of the lowest-energy singlet and triplet excitation of 
2T, 4T and 6T. While for the smaller molecules the first 
exciton stems from the H —>■ L transition, regardless of 
whether we include the exchange and direct e-h interac¬ 
tion, for 6T the lowest-energy excitation, still presenting 
TT-TT* character, is composed by a mixing of H —L (74%) 
and H-1 -A L-l-1 (20%) transitions. The same result is 
obtained also when the BSE exchange term (i7“), and 
hence LFE, are neglected (triplet). From this we can 
assert that e-h correlation effects, given by cru¬ 

cially affect the composition of the first exciton in 6T. 
Fig. [T] also shows the spectrum obtained from an effec¬ 
tive approximation to the e-h interaction term (Hq^^, 
shaded area). This corresponds to considering only the 
G = G' = 0 term of (including a q-independent 

screened e-h interaction). The Hamiltonian then be¬ 
comes -I- ho, where ho is the G = 0 term 

of 
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For details on the derivation of Eq. see Ref. [SH The 
exchange term vanishes, as it contains only the short- 
range part of the bare Coulomb potential v = v — Vq, 
vq being the G = 0 component of v. Excitations re¬ 
sulting from consist of pure single-particle transi¬ 
tions due to which are rigidly red-shifted by the 

“zero-order” e-h interaction from hg (see Eq. .1^ In the 
spectrum of 6T, in addition to the intense peak at 
2 eV, given by H —>• L, a weak peak at about 3 eV stems 
from the H-1 —)■ L-l-1 transition (Table [ll|). The same fea¬ 
ture appears also in the ALDA spectrum (orange line). 
This analogy, supported by the insight from the BSE re¬ 
sults (dominance of confirms that the local and 

frequency-independent cannot properly account 

for correlation effects. In smaller molecules like 2T and 
4T, the OS in the corresponding energy range is concen¬ 
trated only in the first peak, which is indeed largely dom¬ 
inated by H —)• L (Table [TI]) . We conclude that the agree¬ 
ment between BSE and TLDA in the small oligomers is 
determined by the single-particle character of the first 
intense peak. As the size of the molecule increases, cor¬ 
relation effects become more relevant: TDDFT can only 
reproduce the energy of the first excitation but not the 
overall spectral shape. The description of these excita¬ 
tions can be improved to some extent by adopting more 
sophisticated xc functional and kernels, such as global or 
range-separated hybrids. In fact, they improve the un¬ 
derlying density/wave-functions through enhancing the 
exchange term and thus reducing the self-interaction er¬ 
ror. This argument is supported by the findings of Ref. 
1431 where the nature of the lowest-energy excitation of 
6T is analyzed. Within B3LYP, the excitation corre¬ 
sponds to an almost pure H —>■ L transition. Using long- 
ranged corrected PBEO (LC-PBEO), its mixed character 
emerges, with the weight of H-1 —>■ L-l-1 being 7%. This 
kernel leads to a better description of the excitation, by 
reproducing its mixed character. 

We have shown so far that TDDFT can yield, in ex¬ 
cellent agreement with MBPT, the absorption spectra of 
small molecules like 2T, which exhibit a first intense peak 
with mainly single-particle character. In this case, ALDA 
correctly describes the peak and mimics the partial can¬ 
cellation between S — and Eb, providing accurate 
excitation energies. On the contrary, when correlation 
effects play an important role, as in 6T, TDDFT suffers 
from apparent shortcomings. This is, however, not the 
only limiting scenario. The single thiophene ring (IT), 
despite its small size, is another problematic case for 
TDDFT. In Fig. (upper panel), we show the spectrum 
of this molecule computed from ALDA and GoWo+BSE. 
Let us focus on the low-energy part first. It is charac¬ 
terized by two peaks, polarized along the y and x axis 
of the molecule, respectively. Contrary to the case of 
uT (n=2, 4, 6), where the H —L transition yields the 
lowest-energy peak, the first exciton of IT stems from 
a combination of different single-particle contributions, 
including H-1 —>• L (60%) and H —>• L-l-3 (30%). The 
second exciton is cc-polarized and mostly due to H —)■ L 



Energy [eV] 

FIG. 2: (Color online) BSE and ALDA spectra of the IT 
molecule (top), expanded crystal (middle), and crystal (bot¬ 
tom), as shown on the right. Solid lines represent the average 
over the three Cartesian components. In the top panel, dot¬ 
ted (dashed) curves indicate the xx (yy) component of Im 
Em from BSE. LDA and GoWo gaps are indicated by dashed 
lines. A Lorentzian broadening of 0.1 eV is applied to all the 
spectra. 


(84%). This information is summarized in Table HI Our 
BSE results are in g ood agreement with experiments and 
CC calculations,!^ while the ALDA spectrum is poor. 
The first two excitations around 6.5 eV are roughly de¬ 
generate in energy. Through indirect information about 
the single-particle transitions contributing to the ALDA 
excitations, we can attribute this inaccuracy to a lack of 
correlation in TDDFT, in analogy with the previous anal¬ 
ysis of 6T. This is also confirmed in a recent quantum- 
chemical study, showing that only those methods, which 
accurately treat correlation effects, such as CC2, EOM- 
CCSD, and CIS(D), correctly yield energy and charac¬ 
ter of the first two excitations (5.7 and 6.2 - 6.4 eV, 
respectively) 

To better identify the cause of the disagreement be¬ 
tween ALDA and MBPT results, we further analyze the 
spectra in Fig. (upper panel). Compared to the 
spectrum (turquoise shaded area), the onset of triplet ex¬ 
citation energies (green bars) is red-shifted. We recall 
that in both cases the BSE Hamiltonian does not include 



I®* excitation 

2’®'® excitation 


pol. 

composition 

pol. 

composition 

-^0 

X 

H L (100%) 

y 

H-1 L (100%) 

triplet 

X 

H L (99%) 

y 

H-1 L (98%) 

singlet 

y 

H-1 ^ L (60%) 
H L-h3 (30%) 

X 

H ^ L (84%) 


TABLE III: Analysis of the first and second excitation of the 
IT molecule. Polarization direction (pol.) and composition 
in terms of single-particle transitions are provided for 
triplet, and singlet results. 
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FIG. 3: (Color online) Imaginary part of the macroscopic 
dielectric function, averaged over the three Cartesian compo¬ 
nents of IT a) molecule, and b) crystal, computed from BSE 
(top) and TDDFT (bottom). Energies of the lowest-lying 
triplet excitations are indicated by vertical green lines. A 
Lorentzian broadening of 0.1 eV is applied to all the spectra. 
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the exchange term . Hence, and triplet differ 

only by the treatment of the screened e-h interaction: in 
the former case, it is q-independent (see Eq. [^, while 
in the latter one the full dielectric tensor is included.l^ 
Remarkably, the character of lowest-energy excitations is 
the same (Table [TTT| ) . From this we conclude that, no mat¬ 
ter how we treat the attractive e-h term, the lowest-lying 
excitons keep their character. Transferring this knowl¬ 
edge to the singlet spectrum (blue line), we emphasize 
that only the inclusion of the exchange term, yields 
the correct description of the first two excitations. This 
is in contrast to the previous example of 6T, where e- 
h correlation plays the dominant role. In IT, LEE are 
not only responsible for the redistribution of the OS to 
higher energies, but also for the right order and composi¬ 
tion of the lowest excitons. Again, ALDA cannot account 
for these many-body effects, inaccurately describing the 
first two peaks and slightly overestimatin g t he absorption 
onset (see and 1 \TDDFT Table0. As shown in 

Fig. (bottom panel), the main effect of the ALDA ker¬ 
nel (orange line), compared to the independent-particle 
approximation (IPA, shaded gray area), is to blue-shift 
the energy of the first two peaks and to redistribute the 
OS to higher energies. Overall, the TDDFT results re¬ 
flect the discrete spectrum of molecular levels. Therefore, 
above the onset (7 - 10 eV), ALDA captures transitions 
between such localized single-particle states, in better 
agreement with BSE than in the regime dominated by 
strongly bound excitons. 


B. From molecule to solid 

When intermolecular interactions come into play, the 
shortcomings of TDDFT become even more dramatic. 
To extend our analysis in this direction, we go system¬ 
atically to the IT crystalline phase, considered with two 
inequivalent molecules in an orthorhombic unit cell, of 
lattice vectors Oo=9.76 A, bQ=7.2 A, and co=6.67 A.I^ 
We investigate also a model system with lattice vectors 


scaled by a factor 1.25 with respect to oq, bo and cq: we 
refer to this structure as expanded crystal.^ The corre¬ 
sponding spectra are shown in Fig. (bottom and mid¬ 
dle panels). Due to the small size of the molecular con¬ 
stituents, large values of , and hence of F — 14c 

and of Eb, are obtained from MBPT (Table |l]). A sys¬ 
tematic decrease of these quantities, going from the iso¬ 
lated compound to its crystalline phases, indicates en¬ 
hanced screening together with delocalization of the ex- 
citon wave-function, as also found for other molecular 
crystals and polymers.^^HSl A number of intense bound ex¬ 
citons appear in the spectra of both crystal structures. 
These are signatures of the strong e-h interaction in or¬ 
ganic crystals, also beyond the first few excitations. This 
is also confirmed by the spectra shown in Fig. (up¬ 
per panel). By comparing triplet (green bars) and 
excitation energiess (turquoise shaded area), we observe 
a large difference between their absorption onsets, sim¬ 
ilarly to the case of the IT molecule. Also LFE play 
an important role. The singlet spectrum (blue line) is 
blue-shifted by over 2 eV compared to the triplet on¬ 
set and the OS of the low-energy peaks is drastically re¬ 
duced with respect to E[q^^. Given this complexity, it 
is not surprising that TDDFT presents serious problems 
in correctly yielding the spectra. ALDA slightly under¬ 
estimates the absorption onset (Table |T]) and gives only 
two peaks in the respective energy regime, as seen in Fig. 

(bottom panel). Likewise, two intense peaks appear 
also in the IPA spectrum (gray shaded area), ^ 1 eV be¬ 
low the ALDA onset (orange line). In both spectra, the 
continuum starts at about 8 eV, in a region where BSE 
features bound excitonic peaks. This confirms once again 
that ALDA reproduces only “IPA-like” excitations, and 
thus cannot quantitatively capture the spectral features 
of molecular crystals. 


IV. SUMMARY AND CONCLUSIONS 


In summary, through a systematic analysis of the opti¬ 
cal absorption features in oligothiophenes, we have clar¬ 
ified the role of e-h exchange and correlation in describ¬ 
ing optical excitations in molecular materials. For the 
crystal structures, TDDFT turns out to be an inade¬ 
quate approach, being unable to reproduce bound exci¬ 
tons and drastically underestimating the continuum on¬ 
set. In large oligomers, as shown for the case of 6T, 
correlation effects turn up in terms of e-h interaction, 
evidenced by mixed excitations. Conversely, in the sin¬ 
gle thiophene ring, the low-energy excitations are driven 
by LFE. Neither scenario is captured by TDDFT. ALDA 
results are in excellent agreement with BSE merely in the 
case of 2T, where the first intense peak stems from one 
vertical transition. From this we conclude that ALDA 
can be trusted for optical excitations in molecular sys¬ 
tems barely when correlation effects do not play a pre¬ 
dominant role. Only high-level methodologies, such as 
MBPT, which explicitly take into account e-e and e- 
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h interactions, can reproduce optical spectra of organic 
materials, from the gas phase to any of their crystalline 
structure on equal footing. 
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